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We calculate the two-pion wave function for the 7 = 2 S-wave two-pion system with a finite scat- 
tering momentum and estimate the interaction range between two pions. It allows us to examine 
the validity of the necessary condition for the finite-volume method for the scattering phase shift. 
A calculation is carried out with a plaquette gauge action for gluons and a clover-improved Wil- 
son action for quarks at l/a= 1.63 GeV on 32^ x 120 lattice in the quenched approximation. 
We conclude that the necessary condition is satisfied within statistical errors for the lattice size 
L > 32, when the quark mass is in the range m^ = 0.176 — 0.345 GeV^ and the scattering mo- 
mentum in k^ < 0.026 GeV^. We also find that the energy dependence of the interaction range is 
small and it takes 1 .2 — 1 .7 fm for our simulation parameters. We obtain the phase shift from the 
two-pion wave function with a smaller statistical error than that from the conventional analysis 
with the two-pion time correlator. 
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1. Introduction 

The scattering phase shift is an important quantity for understanding a dynamical aspect of 
hadrons. For the 1 = 2 5-wave two-pion system, which is the simplest case, the phase shift has 
been calculated in Refs. ^ ^ The calculations employed the finite-volume method, in which 
the phase shift is related to the energy on a finite volume. It has been proposed by Liischer ^ and 
extended for the non-zero momentum system by Rummukainen and Gottlieb [Q]. 

The derivation of Liischer's and Rummukainen-Gottlieb's formula assumes the condition R < 
L/2 for the two-pion interaction range 7? and the lattice size L, so that the boundary condition does 
not distort the shape of the two-pion interaction. The CP-PACS collaboration calculated the two- 
pion wave function for the ground state of the 1 = 2 5-wave two-pion system and estimated the 
interaction range R from the wave function [Q]. In their case the scattering momentum is highly 
small, ^ ~ 0, thus their work is an examination of the necessary condition for the calculation of the 
scattering length on the lattice. 

In present work we extend their work for the scattering length (the scattering momentum 
~ 0) to that for the scattering phase shift {k ^ 0). For this aim we consider the ground state of 
the system having a non-zero total momentum P = (2n/L)ex in a box satisfying the periodic 
boundary condition. All calculations of this work has been done on VPP5000/80 at the Academic 
Computing and Communications Center of University of Tsukuba. 



2. Finite size formula 

In this section we briefly review the finite-volume method presented by Rummukainen and 
Gottlieb in Ref. [^], with emphasis on the role of the condition for the interaction range. The 
formula has also been derived from another approaches in Ref. ^ and [^. We follow, however, 
from the original derivation in Ref. ||^. 

The two-pion wave function on a finite periodic box of volume is defined by 

<^{x,t) = (0| 7r+(X + x/2,r+f/2) K+{X-x/2,T-t/2) iTrTr;^,?) -e"*'^ -e^'^ , (2.1) 

where \nn;E,P) is an energy eigenstate of the two-pion system with the energy E and the total 
momentum P. 7r+(x,f) is an interpolating operator for 71+ at {x,t). The two exponential factors are 
introduced to remove the trivial exponential factors for the center of mass coordinate X and T. 
In order to relate the wave function to the scattering phase shift, we transform the wave func- 



tion <I>(x,f) in (gj) to that in the center of mass frame <I>cm(x,?) by the Lorentz transformation. 
Here we assume that the two-pion interaction range R is smaller than one half the lattice extent, i.e. 
there exists the region R < \x\ < L/2, where the two pions behave as free particles. In this region 
the wave function satisfies the following two equations [^]. 

{V^ + k^)^CM{x,t)=0 , |-<I>cm(x,0=0, (2.2) 



where k is the scattering momentum related to the invariant mass by ^/s = VE^ — P^ = 2 \J + k'^. 
<I>CM(x,f) also satisfies the boundary condition, 

<J>CM(x,0 = (-l)*^/''''''''"-<J>CM(x + L7[m],0 formeZ^, (2.3) 
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where 7 is the vector operation /[x] = 7X|| +xj^ with the Lorentz boost factor 7 = E j^, X|| 
P(P-x)/p2 and x_L =x-X||. 



e'P- , r = { p I p = (27r/L) • T ' [n] + P/2, n G 



The solution of the ( [2.2| ) under the condition can be given by 

«I>cm(x,0 = -^I -i^ 

per 

'E ' 

per 



(2.4) 



;;73l ^^•Jo(^x) + — •«o('tx)+£Q,„(;t).F/"(il);K^^), (2.5) 
7^ per T' ^tt 



up to overall constant, where i\(kx) is the spherical Bessel and ni{kx) is the spherical Neumann 
function. C/,„ {k) is some constant depending on the scattering momentum k. The first and second 
terms of ( [2.5[ ) consist of the 5-wave component and those coefficients give the 5-wave phase shift 

d{t). 
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tan5o(^) k jL^^^p^-k^' 
This is the Rummukainen-Gottlieb formula [01. 



3. Details of simulation 

In the present work we consider the ground state of the system with the total momentum P = 
and P = {2n/L)tx- We can obtain the scattering length from the energy of the system with P = 



and the phase shift from that with P = {2n/L)tx through the Rummukainen-Gottlieb formula ( P^ ) 
In order to calculate the wave function we consider the correlator. 



F(x,t) = (0| a(x,T) a(P,T,) |0) . 
The operator n(x, t) is defined by 

n(x,T) =££ e'P-^ 7r+(X + ^[x],T) 7r+(X,T) 

R X 



(3.1) 



(3.2) 



The vector operation R represents an element of the cubic group (O/,) for P = and the tetragonal 
group {D^k) for P = {2n/L)ex. The summation over R projects out representation of these 
groups, which equals to the S-wave state ignoring the effects from higher angular momentum Z > 2. 



The operator 11 (P, T^) in (3.1 ) is defined by 



1 



Nr 



a(P,T,) = -^£ [7r+(P,T,;^y)7r+(0,T,;T],-)]^ , 



where 



^e'•P■''^^(x,T,)^;(x) 



75 



(3.3) 



(3.4) 



The operator n^{P,Zs;rij) is defined as ;r+(P,Ts;i^y) by changing ^j{x) to t];(x). The functions 
i^y(x) and Tjy(x) are ?7(1) noise whose property is 



1 Nr 



(3.5) 
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Figure 1: The two-pion wave functions 4>cm(x) at m\ = 0.176 GeV^ for P = and P = {2n/L)tx. The left 
and right panels for each momentum show the wave function on xy- and yz-plane. 



In the present work we take A//; = 2 in (3.3). 



In large T region, we can obtain the wave function for the ground state in (lA) by <I>(x,0) = 
F(x, t) /F{xq, t) introducing the reference position xq. In the present work we set = 20, T — = 
40 and xq = (7,5,2). In the estimation of the interaction range, we analyze the wave function in the 
center of mass frame <I>cM(x,f) at f = 0. It is related to ^{x,t) by <I>cm(7[x],0) = <I>(x,0) outside 



of the interacting region (|x| > R) from the second equation in (2.2). 

Gauge configurations are generated in a quenched approximation with a plaquette gauge action 
at j8 = 5.9 on a 32^ x 120 lattice. The physical quantities are measured every 200 sweeps inde- 
pendently for each quark mass parameter. A clover fermion action with Csw = 1-364 is used. The 
quark propagators are imposed to the Dirichlet boundary condition in the time direction and to the 
periodic boundary condition in the spatial one. The lattice cutoff is estimated as 1 /a = 1.63(5) GeV 
(a = 0.121(3) fm) from the p meson mass. Three quark masses are chosen to give = 0.176, 
0.238 and 0.345 GeV^. The numbers of configurations are 400, 212 and 212 for each quark masses. 

4. Results of interaction range 

In Fig. [T[ we show the two-pion wave functions <I>cm(x) at m| = 0.176 GeV^ for P = and 
for P = (27r/L)ev. The left and right panels for each momentum show the wave function on xy- 
and jz-plane. We find a very clear signal. 

We now consider the two-pion interaction from the ratio, V{x) = V^<I>cm(x)/<I>cm(x). Away 
from the two-pion interaction range, i.e. |x| > R, we expect that V{x) is independent of x and equals 
to —k'^ from (2.2). In Fig. ^ V{x) for the same parameters as for Fig. |l| are plotted. We find a very 
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(a) P = 




(b) P = {2n/L)t, 




Figure 2: y(x) at ml = 0.176 GeV^ for P = and P = [In/Ly^. The left and right panels for each 
momentum show the y(x) onxy- and yz-plane. 



clear signal and V(x) seems to be constant for |x| > 10. We observe a strong repulsive interaction 
at the origin consistent with the negative phase shift of the 1 = 2 two-pion system. 
In order to estimate the interaction range R, we consider, 

U{x) = V2<I>cm(x)/<I>cm(x) , (4.1) 

where is obtained from the two-pion time correlator. According to Ref. we employ the oper- 
ational definition of the interaction range R as the scale where U (x) is sufficiently small compared 
to the statistical error. Strictly speaking, even if |x| takes a large value, U (x) does not vanish and 
has a finite tail. However, the systematic error for the final results of the phase shift due to the 
existence of the tail is buried into the statistical error in this definition. 

In Fig. ||, we show U (x) as a function of |x| for P = and P = (27r/L)ev at the three quark 
masses. We find that the interaction range R takes 

mliG&W^) 0.176 0.238 0.345 

R for P = 13.0 14.0 12.0 

R for P = (27r/L)e:, 10.0 11.0 12.0. (4.2) 

It is at most 14.0 (1.69 fm) and smaller than L/2 = 16. Thus the necessary condition for the 
Rummukainen-Gottlieb formula (2.6) is satisfied within statistical errors at our simulation points. 



5. Results of scattering phase shift 

We estimate the phase shift with the following three methods : 
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Figure 3: U{x) as a function of |x| for P = and P = {2n/L)e.x for several quark masses. 



1. We extract the energy E from the two-pion time correlator and calculate the scattering mo- 
mentum by = {E^ — P^) /4 — m^. The phase shift is calculated by substituting into (2.6). 



The results of the phase shift are shown in the left panel of Fig. ^ (labeled "from T"), where 
the scattering amplitudes, 



A{mji,k) =tan5o{k)/k- \ m\ + k'^ 



(5.1) 



are plotted. 

2. We extract k^ by fitting the wave function <I>cm(x) with the fitting function given in (2.4) 
taking k^ and an overall constant as the fitting parameters. We choose the fitting range |x| > /? 
with 7? given in ( |4.2| ). The results are plotted in the left panel of Fig. |^ (labeled "from W"). 

3. k^ is extracted by fitting V{x) to a constant in the region |x| > /? with R given in (4.2). We 
show the results in the left panel of Fig. ^ (labeled "from V"). 

As shown in the left panel of Fig. ^ the results given by the three method are consistent within the 
statistical errors. The data given by the method ^ ("from V") provides the smallest statistical error, 
so the following analysis is performed with this data. 

In order to obtain the phase shift at the physical quark mass, we extrapolate our results with 
the fit form, 

A{m^,k) =AiQm\+A2Qm\+AQik^+Anm\k^ . (5.2) 
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Figure 4: In the left panel A(m;[,A:) given by the three methods are shown as a function of k^. The dotted 
line is the fit curve for the data given by the method ("from V"). In the right panel our final results of the 



scattering phase shift at the physical quark mass are plotted and compared with the experiment [[10|. 



The fit curves for this fitting are also plotted in the left panel of Fig. 0. As shown in the figure the 
fitting is carried out well. Our final results of the phase shift at the physical quark mass are shown 
in the right panel of Fig.^ and compared with the experiment [[l0|]. Our results are slightly larger 
than the experiment. A possible origin of the discrepancy is finite lattice spacing effects. We must 
leave the confirmation of this to studies in the future. 



6. Conclusion 

In the present work, we have studied the 1 = 2 two-pion wave functions with the scattering 
momentum A;^ ~ and ^ 0. We have estimated the two-pion interaction ranges R from those. It 



has been confirmed that the necessary condition for the Rummukainen-Gottlieb formula (2.6) 
satisfies within statistical errors in our parameters. Moreover, we have estimated the scattering 
phase shift with the two-pion wave function. These methods provide a smaller statistical error 
than that from the conventional analysis with the two-pion time correlator. 
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